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Abstract 

The recent recomputation of the neutrino fluxes from nuclear reactors relaxes the tension between 

q \ the LSND and MiniBooNE anomalies and disappearance data when interpreted in terms of sterile 

neutrino oscillations. The simplest extension of the Standard Model with such fermion singlets is 

the addition of right-handed sterile neutrinos with small Major ana masses. Even when introducing 

(— I \ three right-handed neutrinos, this scenario has less free parameters than the 3+2 scenarios studied 

in the literature. This begs the question whether the best fit regions obtained can be reproduced 



Oh 



by this simplest extension of the Standard Model. In order to address this question, we devise an 

exact parametrization of Standard Model extensions with right-handed neutrinos. Apart from the 
CO 

usual 3x3 neutrino mixing matrix and the 3 masses of the lightest neutrinos, the extra degrees 
CN ' 
^} • of freedom are encoded in another 3x3 unitary matrix and 3 additional mixing angles. The 

parametrization includes all the correlations among masses and mixings and is valid beyond the 

usual seesaw approximation. Through this parametrization we find that the best fit regions for the 

LSND and MiniBooNE anomalies in a 3+2 scenario can indeed be reproduced despite the smaller 

number of degrees of freedom. 
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I. INTRODUCTION 



Despite the success of the standard three-neutrino oscillations in explaining the results of 
solar, atmospheric, reactor and accelerator experiments (see e.g. Ref. |l() by two large and 
one small mixing angles and two distinct mass squared splittings, a number of experiments 
are now hinting to the existence of oscillations at much shorter baselines. This would imply 
the existence of extra mass squared splittings Am 2 ~ 1 eV and therefore the mixing of 
the three Standard Model (SM) active neutrinos with extra sterile neutrino states. While, 
by themselves, none of these hints is very compelling, the fact that they all point towards 
similar regions of the parameter space (Am 2 ~ 1 eV and \U a i\ ~ 0.1) is intriguing. The 
better-known and stronger of these anomalies corresponds to the long-standing 3.8<r excess of 



a. 



v e in the v^ beam observed by the LSND experiment [2j . This anomaly was explored by the 
MiniBooNE collaboration and, while it was not confirmed in the neutrino mode 3j , a similar 
excess, although only at ~ 2a, was found in antineutrinos [4 1 . When interpreted as neutrino 
oscillations, such a difference in the neutrino and antineutrino oscillation probabilities can 
be explained by the CP violation expected to be present in appearance channels. 

The interpretation of these anomalies as mixing between the active SM neutrinos and new 
sterile degrees of freedom 1 is plagued by the tensions between these positive appearance re- 
sults and the constraints derived by the negative neutrino disappearance experiments 6-9] 

nn n 

notably reactor experiments [10h17| for v e and CDHS [18[ and atmospheric 



v^ disappearance. This situation has changec 



neutrino fluxes from nuclear reactors 
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neutrinos for 



recently with the updated computation of the 



201 ] which predict a slightly larger flux, as com- 



pared to the previous estimate, favouring v e disappearance in reactors at 2.2<r and adding 
weight to the sterile neutrino interpretation of the LSND and MiniBooNE anomalies 21]. 



Furthermore, the Gallium neutrino anomaly 
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25| also seems to favour short baseline v e 



disappearance. New fits to a 3+2 neutrino oscillation scheme, i.e., 2 extra sterile neutrinos 
mixing with the 3 SM active ones, show that the tension between appearance and disap- 
pearance data is greatly reduced with the new computation of reactor neutrino fluxes, but 
significant tension remains mainly from the atmospheric and CDHS sector [26j, |27J. Two 



1 Any new light state mixing with the 3 SM active neutrinos needs to be sterile so as not to contribute the 
measurement of the invisible width of the Z |5[, which constrains the number of active neutrinos with 
masses below Mz/2 to 3. 



extra neutrinos are favoured with respect to only one so as to provide the CP violating inter- 
ference necessary to accommodate the negative results in the MiniBooNE neutrino channel 
with the excess found in antineutrinos. More experimental results will be required to clar- 
ify the present tension between the experiments favouring the existence of short baseline 
neutrino oscillations and the null result searches. 

From a theoretical point of view, it is an interesting exercise to consider the particle 
content and interaction Lagrangian that could give rise to the 3+2 scenario that currently 
provides the best (although not-so-good) fit to neutrino data. When these 3+2 mixing 
scenarios are considered in the literature, a general 5x5 mixing matrix is normally assumed. 
Such a matrix can be parametrized by 10 independent mixing angles and 15 phases. Out of 
the 15 phases, 3 can be reabsorbed in charged lepton field redefinitions, 2 are not physical 
and 4 are Majorana phases that could be reabsorbed in the neutrino fields, unless they are 
Majorana particles, and in any case do not play a role in neutrino oscillations. In addition, 
one of the mixing angles is also unphysical, corresponding to an arbitrariness in the definition 
of the sterile flavor states. This leaves us with 9 mixing angles and 6 Dirac phases on top of 
the 5 neutrino masses which are independent and relevant for neutrino oscillations. 

The simplest extension of the SM particle content with singlet fields corresponds to 
the natural addition of right-handed neutrinos. These right-handed neutrinos should be 
singlets of the SM gauge group and allow to write Yukawa couplings between the SM active 
neutrinos and the Higgs field that induce Dirac mass terms after the Higgs develops its 
vacuum expectation value (vev) and thus account for neutrino masses. Moreover, being SM 
singlets, right-handed neutrinos are also allowed to have Majorana mass terms. If we consider 
the addition of 2 right-handed neutrinos, we find that the Yukawa couplings are given by 
a 3 x 2 matrix with 6 independent moduli and the same number of phases. The Majorana 
masses of the right-handed neutrinos only constitutes 2 extra moduli, since a basis in which 
the matrix is diagonal and positive can always be chosen. We therefore see that we have 
only 8 independent moduli and 3 physical phases (one of them Majorana) to be compared 
with the 14 independent moduli (9 angles and 5 masses) and 6 Dirac phases assumed in 3+2 
schemes. Even adding 3 right-handed neutrinos instead of 2, which also seems more natural 
considering the SM particle content, only brings us up to 12 independent moduli (6 masses 
and 6 angles) and 4 Dirac phases. Therefore, the question of whether these simple extensions 
of the SM with a reduced number of degrees of freedom with respect to the general 3+2 



scenario can provide the same best fit for neutrino oscillation experiments naturally arises. 
Notice that, in order to extend the SM by a completely general 5x5 mixing matrix and 5 
independent neutrino masses so as to obtain the 3+2 scenario, a completely general mass 
matrix mixing the 5 neutrino states (3 active and 2 steriles) would be required. This would 
imply, for example, to extend the SM with 3 right-handed neutrinos, so as to allow Yukawa 
couplings and thus Dirac masses plus another extra 4 SM singlets, 2 left-handed and 2 
right-handed, forming 2 Dirac pairs with masses unrelated to the Higgs mechanism and 
thus providing a general 5x5 mass matrix. This particle content seems highly unnatural 
compared to the simple addition of 3 light right-handed Majorana neutrinos. 

Motivated by this question, we have devised a new parametrization of models with extra 



right-handed neutrinos as, for example, the popular type-I seesaw mechanism 28M31] . This 
parametrization takes into account all the relations between masses and mixing parameters 
implied by the reduced number of degrees of freedom discussed above. While several such 
parametrizations already exist in the literature, we find that they have been derived in the 
seesaw limit, i.e., for Majorana masses much larger than the Dirac ones. While this is 
usually an excellent approximation in seesaw models, we found it unsuited for the task at 
hand, since the masses of the new sterile states required for the LSND and MiniBooNE 
anomalies are around 1 eV, not far from the constrained range of the 3 SM neutrinos. Other 



variants of the seesaw mechanism, such as the linear 32J or the inverse 33| seesaws, also 
have very similar Dirac and Majorana masses and thus the standard seesaw approximation 
and the parametrizations derived from it might not be suitable. The new parametrization 
we introduce is exact and makes no assumption on the sizes of the Majorana and Dirac 
mass matrices. In the case in which 3 extra right-handed neutrinos are introduced, the 12 
independent moduli and 6 phases are distributed in a 3 x 3 unitary matrix, which reduces 
to the standard one when the mixing between the sterile and active neutrinos is small, the 
3 masses of the lightest neutrino mass eigenstates, 3 mixing angles which are related to the 
3 masses of the heavier mass eigenstates and an extra 3x3 unitary matrix, making a total 
of 3 masses, 9 angles and 6 phases (4 Dirac and 2 Majorana). 

In Sec. [Til we introduce the new parametrization comparing it with other existing 
parametrizations of the seesaw limit in the literature and discussing the physical ranges 
of the parameters introduced and the relations between them. In Sec. IIHI we apply the 
parametrization to three widely different examples: the seesaw limit, the case of purely 



Dirac neutrinos and the intermediate regime of light Major ana masses that can provide the 
best fit to the LSND and MiniBooNE anomalies, explicitly showing that the 3+2 best fit 
can be reproduced despite the reduced number of degrees of freedom. Finally, in Sec. HVlwe 
summarize the results and draw our conclusions. 

II. PARAMETRIZATION 

Here we discuss a new parametrization of the unknown degrees of freedom of the neutrino 
sector necessary to fully specify the Lagrangian and physical observables when extending 
the SM by fermion singlets (right-handed neutrinos) as in the type-I seesaw mechanism (but 
with arbitrary masses). For simplicity we will here consider the case in which three extra 
right-handed neutrinos are added N^. The Lagrangian would thus read: 

% = JSfsM - ^.(Ma^Ng - (Y N ) ia Nl<pkl + H.c. . (1) 

Here, <p denotes the SM Higgs field, which breaks the electroweak (EW) symmetry after 
acquiring its vev v^w- We have also introduced the Majorana mass allowed for the right- 
handed neutrinos M^ as well as the Yukawa couplings between the neutrinos and the Higgs 
field. The vev of the Higgs will also induce Dirac masses m D = vew^n- Thus, the full 
6x6 mixing matrix L^ot is the unitary matrix that diagonalizes the extended neutrino mass 
matrix: 

UL Utot = , (2) 

\m D M N J yo M J 

where m and M are diagonal matrices. Without loss of generality this diagonalization can 
be performed in two steps: first a block-diagonalization and then two unitary rotations to 
diagonalize the mass matrices of the light and heavy neutrinos, i.e., 



f/«= l ' , (3) 




where U and U' are unitary matrices. We can always choose a basis for the heavy singlets 
such that U' = I (rotations amongst the sterile states are unphysical). The two steps of 
the unitary rotation can be expressed as the exponential of antihermitian matrices, which 
correspond to the Lie algebra of the unitary group. In particular the block-diagonalization 



will correspond to the exponential of a block-off diagonal antihermitian matrix: 
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With a completely general 3x3 matrix. The series of Eq. ((3]) correspond to somewhat 
modified versions of the sine and cosine series, we therefore define: 
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Substituting this expression into Eq. (P2J) we obtain: 



z*U*mU\ 



?*Ms t . 



(6) 



This relation contains all the correlations between the active neutrino masses and mixings 
and the sterile ones and is the starting point of any parametrization. In the seesaw limit 
Mm S> mo so that 



6 ~ m D M^ 



(7) 



represents the generalized 3x3 mixing between the active neutrinos and the heavy mass 
eigenstates and U corresponds to the 3x3 low energy neutrino mixing matrix between the 
flavour states with high accuracy. Then Eq. ([6]) simplifies to the well-known relation: 



U*mU ] = -m D M- l m D . 



(8) 



This last equation has been used to introduce several seesaw parametrizations. For exam- 
ple, the popular Casas-Ibarra parametrization involving an orthogonal complex matrix R 
introduced in 34J exploits the fact that from Eq. (jHJ) the matrix R = %M~ 1 l 2 mBUm~ 1 l 2 has 
to be orthogonal and thus the neutrino Dirac masses can be easily obtained in terms the 
mass eigenstates m and M, the low energy mixing matrix U and the elements of the matrix 
R as mo = —iM l l 2 Rm l l 2 U^. The main advantage of this parametrization is that the heavy 
mass eigenstates in M are part of its set of free parameters. However, the physical range of 
the parameters contained in R can be cumbersome and lead to complications in the scan of 
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the seesaw parameter space. Indeed, R can be parametrized by three complex angles with 



unconstrained imaginary parts, see Ref. 35[ for a detailed discussion 



Another convenient and popular parametrization, discussed in detail for example in 
is obtained by diagonalizing the Dirac mass matrix through two unitary rotations: tjid = 
UrDU l where D is a real diagonal matrix and Ur and Ul are unitary. Substituting in 
Eq. ©: 

D^UI^ttiU'UlD- 1 = -U^M^Ur. (9) 

Thus, the heavy mass eigenstates M and the mixing Ur can be parametrized in terms of 
the light masses and mixings m and U together with a unitary matrix Ul and the three real 
numbers contained in D by performing the diagonalization of Eq. Q. This parametrization 
replaces the complex orthonormal matrix of the Casas-Ibarra parametrization, R, by a 
unitary matrix, Ul, allowing more convenient scans of the parameter space. Conversely, the 
heavy mass eigenstates are no longer free parameters but derived. 

A generalization of this second parametrization away from the seesaw in which Mn ^> ttld 



limit was performed in Ref. 37]]. In that case, the Lagrangian parameters, i.e., D, Ur, Ul 
and Mjv where used as free parameters. While this choice allows a very simple reconstruction 
of the Lagrangian, relating these parameters with the mass eigenvalues and mixing angles 
can be cumbersome and we therefore propose a different approach. 

A. Parametrization outside the seesaw limit 

In order to make an exact parametrization, we start with Eq. fl6]) and introduce a biunitary 
diagonalization of 0: 

e = v L ovl (io) 

where 9 is a positive diagonal matrix of angles. Under this diagonalization s = Vl sm(9)V^, 
c = VlCOs(9)VI and c = VrCOs(9)Vr. Thus, Eq. (jHJ) becomes: 

cot^V^fTmf/ViCot^) = -V&MVr. (11) 

Therefore, the masses M and the mixing Vr can be obtained in terms of the masses m and 
mixings U, together with a unitary matrix Vl and the three real angles contained in 9 by 



performing a diagonalization. The full mixing matrix is then given by: 



£/, 



lot 



V L cos(6)VlU V L sm(e)Vk 
-V R sm(e)VlU V R coa{9)v£ 



(12) 



Regarding the independent phases contained in the matrices Vl and U, parameter counting 
tells us that the original Lagrangian had 9 independent phases, since all the elements of 
the Dirac mass matrix mr> can be complex. Of these 9 phases, 3 can be absorbed by field 
redefinitions of the charged leptons. This is not an option for the neutrino fields since they 
have Majorana mass terms. To identify the remaining 6 independent phases we write down 
the three unitary matrices in the form <l> 1 V r $2 5 where $i and $2 are diagonal matrices of 
phases and V has the usual CKM form. 
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y -si 3 e td C13 J 

Notice that, from the definition of Vl and Vr in Eq. fllOp . the phases contained in $ 2 for 
the matrix Vl can be absorbed in the definition of Vr. As for the matrix U, the phases 
contained in $2 are the usual Majorana phases of the 3x3 neutrino mixing matrix, while 
those contained in $1 either appear in combination with the $1 of Vl or can be absorbed 
via the field redefinitions of the charged leptons. Therefore, only the $1 of either U or 
Vl is independent. We then choose the following parametrization for the two independent 



unitary matrices U = V(9 i2 , #23, 6> 13 , 5) ■ diag(l, e~ ia2 , e~ ta3 ) and V L = diag(l, e~ ia ^ , e" m s ) • 
V(9i2, $23) ^135 ^ L )- Notice that a 2 and a 3 can be considered Majorana phases, since they 
would become unphysical if the active neutrinos were Dirac fields, while 5, 5 L , a\ and a 3 
correspond to Dirac phases. 

An alternative approach to the parametrization suggested above is to perform the equiv- 
alent of the Casas-Ibarra approach outside of the seesaw limit. In this case, we introduce 
the matrix 

B = iy/M-^s^cUy/m, (14) 

which by virtue of Eq. ([6]) is complex orthogonal. From this definition we can do the 
computation 

F = tU^B T VM~ 1 = c _1 s = V L tan(0) v£, (15) 
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where U, M, m and B are choosen as the independent set of parameters, while Vl, Vr and 
are quantities that can be computed from the bidiagonalization of F. Once all the quantities 
are known, the full mixing matrix C/totj and therefore also the fundamental parameters in 
the Lagrangian, can be computed using Eq. ( TT2J) . The advantage of this parametrization 
is that it contains all of the physical masses as parameters, while the disavantage is that it 
does not have a clear interpretation of the complex orthogonal matrix B (while in our first 
parametrization, is directly related to the active-sterile mixing). 

In what follows we will use the first parametrization, whose independent parameters 
consists of m, U, Vl, and 0. 

B. Physical ranges 

We will now discuss the physical ranges that the angles and phases that define the neutrino 
mixing matrix should take in order to reproduce all the distinct physical situations avoiding 
overcounting. Notice that, apart from the freedom that allows us to redefine the charged 
lepton fields up to an overall phase used to absorb three of the phases, in the neutrino 
fields we can also perform the transformation v± —¥ —Vi and their Majorana masses remain 



unchanged. Therefore, following Ref. 38[, two sets of mixing parameters Oij, 5, «j, Of ,-, 5 L , 
af and 0\^ 5', a\, 0'^, 5 /L , a'f will describe the same physical situation if the corresponding 
mixing matrices are related through the transformations: 

U tot (0[,p *', aj, 0%, 6' L , af) = P'UtotiOij, 5, a u Ofj, S L , af)P n , (16) 

where P l and P n are diagonal matrices with either -1 or 1 in the diagonal and represent the 
corresponding field redefinitions of charged leptons and neutrinos respectively. Moreover, 
only the first three rows of U to t i n Eq. (TT2|) are physical, since only the active neutrinos take 
part in charged current interactions. Therefore, Eq. (fTBl reduces to: 

V[cos(0)V' L ] U' = P l V L Pcos(0)PV£P l P l UP n 

V' L sin(#)V^ = P l V L P sm(0)PV^P s , (17) 

where, for simplicity, we have omitted the explicit dependence of the matrices on the angles 
and phases defining V[ = VhiO'hi 5 /L , at' i L ) and so on. P and P s are arbitrary matrices with 
the same structure as P l and P n . Thus, transformations of the form 

V[ = P l V L P 

9 
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TABLE I: Set of transformations of angles and phases that describe the same physical configuration. 
The restrictions that each transformation allows to set on the phases and mixing angles in order 
to bring them to their physical range is also displayed. 



V' = P S V R P 



R 



U' = P l UP n , 



(18) 



represent the same physical situation. Since Vr is not independent from Vl and U it is 
important to confirm the consistency of the set of transformations of Eq. ffT8|) with Eq. (fTTj) . 
Following Ref. 38j we show that with different choices of P , P n and P a set of transfor- 
mations among phases and angles can be denned such that the physical range of the 6 angles 
can be chosen to be [0,7r/2], of the 2 Majorana phases (— 7r/2,7r/2] and (— 7t,tt] for the 4 
Dirac phases. Tabled] shows the list of transformations among the angles and phases that can 
be used to bring them to their physical ranges and the corresponding values that P l , P n and 
P take for each of them. Notice that when P l or P n transform, the corresponding charged 
lepton or neutrino field redefinitions are required to recover the original configuration. 
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III. APPLIED EXAMPLES 

In order to get a more intuitive handle on the proposed parametrization, we will now 
examine how it applies to different known models. In particular, we will discuss the limiting 
cases of the type-I seesaw limit and Dirac neutrinos, as well as the intermediate regime of 
light Majorana masses that can provide the best fit to the LSND and MiniBooNE anomalies, 
explicitly showing that the 3+2 best fit can be reproduced despite the reduced number of 
degrees of freedom. 

A. The seesaw limit 

In the seesaw limit Mn ^> rrio we can expand G to first order and obtain Eq. ([7j). In 
addition, the basis where U' = I corresponds to the basis where M^ is already diagonal also 
to first order, such that ~ rripM -1 . Furthermore, since is small, we have c ~ 1 and 
s ~ 0. Using this last relation, it is easy to deduce 

V L sin 2 (0)y£ ~ V L 6 2 Vl = 00 f = m} D M- 2 m D . (19) 

The last equality corresponds to the coefficient of the d = 6 operator that is obtained in the 



type-I seesaw after integrating out the heavy neutrino degrees of freedom [39J . Thus, in this 
limit, the extra degrees of freedom 9 and Vl are simply obtained from the diagonalization 
of the d = 6 operator. The remaining parameters, i.e., m and U are contained, as expected 
in the coefficient of the Weinberg d = 5 operator from Eq. (jSJ), responsible for the light 
neutrino masses. Thus, we find that, as expected, the d = 5 and d = 6 operator suffice to 



reproduce a 



matrices |39| 



1 the parameters of the original Lagrangian, as well as the full mass and mixing 



Comparing with the Ul parametrization of Ref. 36J, D = U^thdUl, we have 



M- l m D = M~ l U R DUl = f ~ V R 9Vl (20) 

Thus, the 8 of our parametrization is intimately tied to D through this relation, and so are Vl 
and Vr with Ul and Ur, with the difference that, in this Seesaw limit, our parametrization 
corresponds to bidiagonalizing M~ l mD rather than m^. It is also easy to relate 9 to the 
Casas-Ibarra parametrization through the R matrix in a similar fashion, expressing mo, and 
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thus 0, as = —iM 1 l 2 Rrn L l 2 XJ 1[ . The procedure to o 
operator, i.e. from 9 and Vl, can be found in Refs. 40 



Dtain the matrix R from the d = 6 



4l|- 



B. Dirac neutrinos 



In the case when M^ = 0, we recover the limit of Dirac neutrinos. In this limit, the full 
mass matrix has the form 

M D = 



mjy 



m D 



(21) 



and can be diagonalized as 



Ul t M D U tot 




(22) 



where 



Utot ~ 7! 



v L D vE 



(23) 



and m = VE T i r n^ ) VE is diagonal. However, VE is an arbitrary rotation of the sterile right- 
handed component of the Dirac neutrinos and we can therefore choose a basis where VE = 1- 
In the parametrization we propose, this is reproduced in the case where 9 = ix /4, which 
corresponds to maximal mixing between the laft and righ-handed fields, as expected. In this 
case, Eq. ( TTTT) turns into 



VTU*mU ] V L = -VZMV, 



R 



■R, 



(24) 



with the obvious solution Vr = WVl and M = —m. Inserting this into Eq. ffTUj) for U tot , 
we obtain 



t/tot = -^ 



U 



rt 



(25) 



R I - —I U U 

! I -y R K f f/ 1 )' ' V2 ^_i l 

which exactly corresponds to the scenario of Dirac neutrinos. It should be noted that there is 
an arbitrariness in U = VlV^, since any combination of Vl and Vr satisfying this condition is 
equally valid. This is related to the fact that the rotations VE of the right-handed neutrino 
fields are unphysical and thus only the physical mixing matrix U remains. 
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0.160 
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1.64 


0.47 


0.88 


Ref. [27] 


0.130 


0.130 


0.134 


0.080 


1.52 


0.9 


1.6 


NH2 


0.128 


0.135 


0.134 


0.080 


1.52 


0.9 


1.6 


IH2 


0.130 


0.130 


0.137 


0.080 


1.52 


0.9 


1.6 



TABLE II: Best fit points for the 3+2 scenario to short baseline neutrino oscillation data performed 
in Refs. 



26 



27l | compared with some sample points chosen to reproduce these values with either 
a normal or inverted hierarchy pattern. The values of the parameters for these sample points are 
shown in Tab. Mil Am^ and Am^ are given in eV 2 . 



C. LSND/MiniBooNE 3+2 best fit 

We will now address the question of whether the best fit found in the 3+2 scheme with 
a general 5x5 mixing matrix for the short baseline neutrino oscillation experiments can be 
reproduced with the addition of only 3 right-handed neutrinos to the SM particle content. 
Notice that, while the general 3+2 scheme has a 5 x 5 mixing matrix with 9 independent 
mixing angles, 6 Dirac phases and 5 masses that can play a role in neutrino oscillations, the 
only relevant parameters are the 4 mass square splittings, the elements of the 3x3 submatrix 
relevant for solar and atmospheric neutrino oscillations, and 4 moduli and a phase for the 
short baseline oscillations that could accommodate LSND and MiniBooNE. Two studies 
have been performed so far in the 3+2 framework including the new reactor fluxes, obtaining 
slightly different results for the best fit values of the extra mass splitting and elements of the 
mixing matrix [26|, |27|. The respective best fits are reported in Tab. [Til Notice that, apart 
from the two extra splittings Am^ and Am^, only the moduli |£/ e4 |, {U^l, |f/ e5 | and {U^l 
are relevant for short baseline neutrino oscillations together with the CP violating phase 
<p = arg([/ e 4[/* 4 [/* 5 [/ M 5) that can accommodate the null results of MiniBooNE in neutrinos 
and the positive signal in antineutrinos. 

The possibility of explaining the LSND anomaly through the inclusion of light right- 
handed Majorana neutrinos was discussed in Ref. 42|, while the first steps towards answering 
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the question of whether the best fit of Ref . [26J could be reproduced by adding 2 right-handed 
neutrinos to the SM Lagrangian were already given in Ref. 37j. There it was shown that, 
when assuming all mixing elements real, the values of \U e 4\, \U^\, \U e5 \ and \U^\ were in 
the correct ballpark compared to those required for short baseline neutrino oscillations when 
obtaining the correct pattern of masses and mixings in the solar and atmospherics sectors 
with an inverted hierarchy pattern. Therefore, it is expected that, when extending the 
analysis to the complex case introducing phases, the best fit would be possible to obtain. 

Here we will extend the analysis to the addition of 3 right handed neutrinos and show 
that both the best fits from Refs. |26l . 1271 ] can be obtained through our parametrization while 
at the same time reproducing the correct solar and atmospheric masses and mixings pattern 
with either a normal or an inverted hierarchy. We have performed a scan of the parameter 
space thr oug h Markov Chain Monte Carlo techniques (MCMC) through the MonteCUBES 



software 



43j and selected 4 points that reproduce well the best fits from Refs. [26|, |27| for 



either a normal or an inverted hierarchy. The values for the relevant parameters for short 
baseline neutrino oscillations are shown in Tab. [Ill and compared to the best fit values from 
Refs. [26|, |27j. The points are labelled either with NH or IH depending on the hierarchy 
assumed and 1 or 2 when trying to reproduce the best fit from Ref. 26J or Ref. 27] respec- 
tively. The elements of the U matrix (except for the phases) were fixed to the current best 
fit points of solar and atmospheric neutrino oscillations from Ref (l|, since we found that 
3x3 mixing matrix that governs these oscillations, i.e. VI cos(9)VlU, did not deviate much 
from U and we always obtained its elements within the la region. The Am^ and Am^ 
splitting were similarly fixed to their best fit values assuming either a normal or an inverted 
hierarchy. We also fixed #13 = 0.1, a relative large value motivated by the recent hint of 
sizable #13 from the T2K experiment 44j. The rest of the parameters were kept free in the 
MCMC and their values are shown in Tab. IIIII for each of the 4 points chosen. 

The value of the third mass splitting Am^ is controlled by the value of the angle 83. 
We chose a small value in all cases that tends to make this third neutrino heavy and its 
mixings small, so that it mostly decouples. In the examples displayed it was close to the 
keV scale. However, we found that changing this parameter does not affect much the rest 
of the values. Thus, decreasing # 3 , the third neutrino can be made much heavier and make 
it decouple while still keeping good agreement with the best fits of Refs. 26|, |27J]. One 
possibility is to choose #3 ~ 0.1° — 1°. In this case, the third neutrino will have a mass 



14 





NH1 


IH1 


NH2 


IH2 


mi (eV), 5, 5 L 


0.02, 38°, 126° 


0.01, -132°, 132° 


0.02, 12°, -168° 


0.02, 160°, 166° 


«2, «3 


-16°, 42° 


-88°, -4.9° 


39°, 34° 


84°, 38° 


nL aL qL 
a 12i ^23' p 13 


78°, 52°, 45° 


58°, 47°, 56° 


61°, 27°, 27° 


64°, 33°, 55° 


«2 . a 3 


-120°, -118° 


-45°, 147° 


-77°, -74° 


-46°, -82° 


#li #2, #3 


9.8°, 18°, 1.6° 


12°, 17°, 1.3° 


6.2°, 13°, 3.0° 


11°, 13°, 0.31° 



TABLE III: Values of the free parameters in the MCMC scan corresponding to the 4 sample points 
chosen. The mixing angles of the U matrix as well as m<i and 7713 were fixed to their present best 
fit value for solar and atmospheric neutrino oscillations taken from Ref. [l|]. 



around a few keV and mixings of order 10~ 4 — 10 -3 with the active neutrinos, close to the 
region in which dark matter could be explained by sterile neutrinos (see, e.g., Ref. [45]). 
On the other hand, decreasing #3 to smaller values increases the mass of the sixth neutrino 
as l/(# 3 ) 2 and reduces its couplings to the active ones as #3, making it decouple. Notice, 
however, that its contribution to light neutrino masses, i.e. to the Weimberg d = 5 operator, 
does not decouple in this manner, since the combination m a b = U* 6 MqU^ of Eq. (jSJ) remains 
constant when decreasing #3. Therefore, the analysis performed here has extra degrees of 
freedom with respect to the adition of only 2 right-handed neutrinos considered in Ref. 37J 
and, hence, obtaining the best fit values of Refs. 26|, |27J is easier. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have addressed the question of whether the best fit points found for 
short baseline neutrino oscillations trying to accommodate the MiniBooNE and LSND re- 



sults of Refs. 



26 



271 ] can be reproduced with the more economical addition of 3 Majorana 



right-handed neutrinos to the SM. Such a scenario not only looks more natural from the 
point of view of the particle content, making the lepton and quark sectors symmetric and 
accommodating neutrino masses with a minimal extension, but also involves less degrees 
of freedom than the extension of the 3x3 neutrino mixing to 5 x 5. Indeed, the general 
5x5 unitary mixing matrix of the 3+2 schemes contains 9 independent angles, 6 Dirac and 
4 Majorana phases plus the 5 distinct neutrino masses, while the addition of 3 Majorana 
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right-handed neutrinos only involves 9 independent angles, 4 Dirac and 2 Majorana phases 
and 3 masses (the other 3 being dependent combinations of the other masses and mixings). 

In order to address this question we have developed a new parametrization of neutrino 
masses and mixings when the SM is extended by Majorana right-handed neutrinos. The 
advantage of this parametrization is that it is exact and, thus, valid for all regions of the 
parameter space. In particular, away from the seesaw limit that most existing parametriza- 
tions assume. This parametrization contains all the relations among masses and mixings 
defining a set of independent parameters from which all others can be easily derived. It can 
therefore be very useful for any scenario in which the seesaw relation, i.e., M^ ^> mp, is 
not met with enough precision to justify the standard approximations as, for example, some 
cases of inverted or linear seesaw models. The parametrization is of course also valid for the 
standard seesaw and its correspondence to other popular parametrizations for this limit has 
been discussed. 

For the purpose at hand, the new parametrization allowed us to scan the parameter 
space down to very low Majorana masses, searching for points close to the best fit values of 
Refs. [26|, |27| . Four such points were chosen from a MCMC scan as samples to show that 
indeed the addition of 3 Majorana right-handed neutrinos can reproduce the more general 
3+2 fits. These four points correspond to the two best fits available in the literature and the 
two mass hierarchies, normal and inverted, which currently can fit solar and atmospheric 
neutrino oscillations. All four samples are very close to the corresponding best fit regions 
showing that they can be accommodated for any choice of the mass hierarchy. Moreover, the 
mass of the third extra neutrino mass eigenstate, not necessary for the fit of short baseline 
neutrino data, can be chosen at any scale and still good fits are obtained. Thus, it could be 
used to try to account for the observed dark matter component of the Universe as a sterile 
neutrino candidate or made very heavy. 
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